The low-energy excitation spectrum of one-dimensional dipolar quantum gases 
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We determine the excitation spectrum of a bosonic dipolar quantum gas in a one-dimensional 
geometry, from the dynamical density-density correlation functions simulated by means of Reptation 
Quantum Monte Carlo techniques. The excitation energy is always vanishing at the first vector of the 
reciprocal lattice in the whole crossover from the liquid-like at low density to the quasi-ordered state 
at high density, demonstrating the absence of a roton minimum. Gaps at higher reciprocal lattice 
vectors are seen to progressively close with increasing density, while the quantum state evolves into 
a quasi-periodic structure. The simulational data together with the uncertainty-principle inequality 
also provide a rigorous proof of the absence of long-range order in such a super-strongly correlated 
system. Our conclusions confirm that the dipolar gas is in a Luttinger-liquid state, significantly 
affected by the dynamical correlations. The connection with ongoing experiments is also discussed. 



Introduction.- Ultracold quantum gases with dipolar 
interactions are currently being produced in laboratory, 
where atomic 52 Cr atoms have been Bose-condensed [lj, 
2j, following earlier theoretical predictions [HQ. Exper- 
iments have been suggested [|[ aimed to produce molec- 
ular gases with large dipolar strengths, and a few labo- 
ratories worldwide are working along these lines. In fact, 
dipolar quantum gases are emerging as competitive real- 
izations of quantum devices Q and as a laboratory for 
investigating strongly correlated regimes @, [U and novel 
quantum phases [9|, [lfj , in which quantum fluctuations 
are enhanced by exploiting techniques acquired for an ac- 
curate manipulation of atomic gases. These include the 
possibility of lowering the temperature; of tuning the in- 
teractions in both their long-range tail [r"J ] and in the 
strength of their short-ran ge par t [l2T ] by means of the 
Fano-Feshbach mechanism [lj, [l4[ to let emerge the dipo- 
lar character [l5[ ; of reducing the dimensionality down to 
one (ID), as already performed in other systems [H, H3|- 

One-dimensional quantum gases are naturally inclined 
to be strongly correlated [l8[ [2^] . We have more recently 
predicted that dipolar bosonic quantum gases confined in 
quasi- ID geometries can reach correlation regimes well 
beyond those of the (already strongly correlated) Tonks- 
Girardeau (TG) gas crossing over to a Dipolar- 

Density- Wave (DDW) state at very large densities n on 
the scale of the potential range, where the atoms ar- 
range into an ordered state regularly spaced by n _1 [3fj| . 
By a back-to-back comparison with Reptation Quantum 
Monte Carlo (RQMC) [2(| simulational data, we have 
shown that at the level of the static structure factor 
the crossover can be described by a Luttinger-liquid the- 
ory with exponent K < 1 continuously decreasing from 
K = 1 nro — > to K — > as nro — > oo. Finally, we have 
predicted the corresponding signatures in the collective 
excitations of the trapped gas [211 ] . 

Beyond the evidence emerging from the static struc- 



ture of the fluid, a clear-cut demonstration of Luttingcr 
behavior requires further understanding of the excita- 
tions in the homogeneous dipolar gas. In particular, 
answers to two relevant questions are not obvious from 
the beginning. First, whether roton-like excitations may 
show up in the dipolar gas at finite wavevectors. Sec- 
ond, whether the quantum fluctuations of the phonon 
field prevent the existence of long-range order at large 
densities, namely whether the crystal order parameter 
vanishes in the thermodynamic limit. In fact, exploit- 
ing the uncertainty-principle instead of the Bogolubov 
inequality, Pitaevskii and Stringari [22| have worked out 
an extension of the Hohenberg-Mermin- Wagner theo- 
rem [23l [2^ which yields more accurate upper bounds 
to the size of the order parameter at zero temperature, 
where the quantum fluctations dominate. When applied 
to specific systems, the inequality may allow to rule out 
the existence of long-range order, as in the case of, e.g., 
ID antiferromagnets and crystals (22j. Both questions 
above would have a definite answer if the system were in 
a Luttinger-liquid state, for which there is no long-range 
order nor roton minimum. 

In this Brief Report, we find that this is indeed the 
case, after computing by RQMC the low-energy exci- 
tation spectrum up to eight reciprocal lattice vectors 
G m /n = 2nm in the whole crossover. The evolution of 
quasi-long-range order from the TG to the DDW state 
emerges as a progressive closing of the gaps in the excita- 
tion spectrum with increasing the order m. By the same 
token, we demonstrate the absence of a roton minimum 
at 2ir in the whole crossover and that dynamical effects 
play a significant role in building the Luttinger state. Our 
results, analyzed by means of the uncertainty-principle 
inequality [22J , also rule out the existence of long-range 
order in this super-strongly correlated quantum gas and 
confirm that the ID dipolar gas is in a Luttinger-liquid 
state. 
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The model and the RQMC method.- We model the ID 
dipolar Bose gas by considering TV atoms with mass M 
and permanent dipoles moments arranged along and or- 
thogonal to a line, yielding purely repulsive interactions. 
The Hamiltonian is 

H= —V— 1 \ - 1 (l) 

r? 2— 1 dx 2 r 3 . 2— 1 \xi - xA 3 ' 

* i 1 " i<j ' J 1 

in effective Rydberg units Ry* = % /(2Mr%). The effec- 
tive Bohr radius r = MCdd/ (2irh 2 ) is expressed in terms 
of the interaction strength Cdd = VaP^ for magnetic and 
Cdd = d 2 /eo for electric dipoles [31]]. The dimension- 
less parameter r s = l/(nro) determines the interacting 
regime in terms of tq and of the linear density n. Since 
the potential-to-kinetic energy ratio scales as l/r s = nr , 
large densities yield to strong correlations, at variance 
with Coulomb systems. 

We determine the ground-state properties and the exci- 
tation spectrum by resorting to the Rcptation Quantum 
Monte Carlo technique [20]. This is in essence a path- 
integral method at zero temperature, where the ground- 
state distribution is directly sampled in the internal part 
of the path. Thus, the computation of the structure of 
the fluid and of the imaginary-time correlation functions 
for suitable long projection times is conceptually straight- 
forward and practically easy, since possible biases arising 
from mixed averages are ruled out by definition. In par- 
ticular, from an analysis of the imaginary-time density- 
density correlation function we determine the low-energy 
excitation spectrum while the parameter nro spans the 
whole crossover from the TG to the DDW state. 

We use a trial wave-function that is a product of two- 
body Jastrow factors iptriai — Yii<j e u ^ Zi ~ z A\ As we are 
interested in long-range behavior, we actually take the 
Luttinger-liquid expression 

Atrial OR) OC JJ | Sin j- (Xi - Xj ) | l/K . (2) 

i<j 

which in the low-density limit implies K = 1 @ and 
recovers the wave-function of spinless non-interacting 
fermions. Different choices of the wave functions, such 
as the product of gaussians centered on the lattice sites 
R m = mn , result into different time-step extrapola- 
tions, but eventually lead to negligible differences in the 
computation of the static and dynamic structure factors. 

A few technical details are in order. We perfom simu- 
lations for different values of the number N of bosons in 
a square box with periodic boundary conditions, namely 
N = 40, 60, 80, 100, reaching in selected cases N = 200. 
We check that we are able to take care of finite-size effects 
by summing the interactions over ten simulation boxes. 
Finally, the energies are extrapolated to their thermody- 
namic limit after removing the time-step dependence. 

The resulting energy per particle £(nr n ) / Ry* as a func- 
tion of nro has been provided in [2l|, [25[ , together with 
an accurate analytical form of it, useful for further ap- 
plications. We remark here that e(nro)/ Ry* recovers the 





> 1 




\ % A 

a — J\ j \ j 







2n 4% 671 8tc 
q/n 

FIG. 1: S(q) in dimensionless units for a dipolar gas with N = 
40 particles and different values of nro = 0.01, 50 and 1000. 
Decreasing slopes as q — > and the emergence of additional 
peaks correspond increasing nro values. 

known limiting behaviors e(nr ) / Ry* ~ (ir 2 /3)(nr ) 2 for 
nr <C 1 in the TG regime and e(nr )/Ry* ~ C(3)(w ) 3 
for nro 3> 1 in the DDW limit We set from now on 
units of n _1 for lengths and of Ry* for energies. 

Low energy excitations from dynamical structure 
factor.- We determine the low-energy excitations after 
computing the imaginary-time correlation function of the 
density operator [26[ p q = exp(—iq- r;(r)) creating a 
density fluctuation with wave vector g, that is F(q, r) = 
(p q (r)pJ(0))/iV, where the sum in p q spans over the 
number of particles N located at position fl and r is 
the imaginary time. F(q,r) is related to the dynamical 
structure factor by F(q,r) = J Q dtu exp(—u!T)S(q,uj), 
yielding the static structure factor for t — 0, namely 
S(q) =F(q,0) = f™cLjS(q,cj). 

To give a system overview, we first display the S(q) in 
Fig. Q] as previously determined Q by RQMC for vari- 
ous densities between the TG and DDW limits. Peaks 
at q/n = 2irm (m integer) progressively disappear with 
decreasing density as the dipolar gas approaches the spin- 
less fermionic liquid. 

We track this smooth evolution by investigating the 
low-energy excitations as extracted from the dynamical 
structure factor S(q,oj). From the expression 

S(q,u,) = J2\(n\p q \0)\ 2 ^~"n), (3) 

n 

we estimate the energy dispersion of the collective exci- 
tations by fitting the imaginary-dependence of F(q, t) as 
a sum of exponentials J ^ li Ai(q)e UJi ^ T corresponding to 
multiple modes. Then, the fit yielding the best \ 2 value 
is chosen. 

Fig. [^displays the resulting RQMC energy dispersion 
u)(q), for the case with N = 40 at nr = 1, 10 and 1000 
namely in the low, intermediate and very high density 
regimes. In spite of the finite size effects, the overall 
qualitative behavior is already clear. The phonon soft- 
ens at low g-values, while the density decreases. On the 
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FIG. 2: Lowest excitation energies u){q) in Ry* units and 
scaled by (nro) 2 , for a dipolar gas with N — 40 and different 
values of nro = 1, 10 and 1000 as in the legend. Symbols with 
error bars represent the RQMC data extracted from (5). the 
solid line is a guide to the eye. The curve at nro = 1000 is de- 
pressed by a factor of 5 for graphical reasons. Inset: zoom on 
the uj(q) at nr — 1000 up to q/n = 87r for different F(q,r) 
models: multimode (solid) and Feynman (dashed) approxi- 
mation. Dotted line: periodic replica of the first bump. 



other hand, the gap at 2ir seems to be always closed at 
all densities. As represented in the inset by the dotted 
curve, even at nro = 1000 the RQMC excitation spec- 
trum (solid line) is very different from what could be 
obtained by replicating the portion from q = to q = 2n 
(dotted line), indicating a non periodic structure. The 
dashed line instead, represents to(q) obtained from the 
Feynman relation Lu(q) = e(q)/S(q), which provides only 
an upper bound to uo(q) in terms of the static structure 
factor S(q) and of the kinetic energy e(q) = H 2 q 2 /2m. 
Indeed, the Feynman relation is clearly seen to overes- 
timate the RQMC value for q > 2tt, while it seems to 
account for the io(q) between q = and q = 2ir. Con- 
sidering that the Feynman relation is expected to yield 
better results as q/n — ► 0, we can also anticipate that at 
lower densities the q range where the Feynman relation 
is reliable will shrink (see below). 

A quantitatively reliable measure of the gap sizes re- 
quires an accurate size effect analysis. Fig. [3] displays 
the 1/N scaling of w(q/n = 2tt) for nr = 0.01, 0.1, 1, 10 
and 1000. The fit to the RQMC data (symbols in the 
figure) yields the linear scaling u>jv(<7 = 2n) = c(nr )/N 
with the constant c(nro) being an increasing function of 
nro. Thus, u{q/n = 2ir) — > as 1/N — > 0, the gap is 
clearly closed at all densities, demonstrating the absence 
of a roton minimum. This is not the case if we use the 
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FIG. 3: 1/N scaling of ui(q — 2tt) in arbitrary units at nro = 
0.01, 0.1, 1, 10 and 1000 as in the legend. Symbols with error 
bars: RQMC data. Solid lines: fit to the data. Inset: zoom 
of the RQMC data and related fits for the cases with nro = 1 
and 1000. Filled symbols and dashed lines: Feynman relation 
(single-mode approximation). Open symbols and solid lines: 
two- mode analysis yielding the best \ 2 - 



Feynman relation, corresponding to a single-mode ap- 
proximation in Eq. ([3]) [201 ] . At the intermediate-to- low 
density nro = 1 the fit yields a finite gap value (see the in- 
set). Since the Feynman relation is built up from static 
quantities, we may conclude that dynamical effects, as 
embodied in the multimode analysis of ([3]) play a signif- 
icant role, leading to qualitatively different conclusions. 
A similar multimode analysis performed at q = 27rm, 
m > 2 shows the existence of open gaps, which progres- 
sively close while the quasi-ordered state is approached. 

Absence of long-range order.- Using these results, we 
can derive a strict upper bound for the order parameter 
of the solid pq = N~ 1 (^2 m exp(iG • r m )) with G a vector 
of the reciprocal lattice, and rigorously test the quali- 
tative conclusions from the inset of Fig. [21 namely that 
no long-range order may exist in our ID dipolar quan- 
tum gas. We closely follow the derivation of Pitaevskii 
and Stringari [22| . By applying the uncertainty-principle 
inequality {A\A){B\B) > \([A^B)}\ 2 to the oper- 
ators A = px + Q and B = dp^/dt one has S(q + 

G)Jdcju; 2 S(q,Lj) > j^p 2 d (q- (q + G) 2 )). From the 

RQMC data we know that as q -v 0, S{q + G) -> |g| 2 ^ -1 
, while the second moment of S(q,ui) vanishes as |g| 3 . 
Thus, the order parameter in the long wavelength limit 
vanishes as p% < q mm i 2K ^) with K > 0. Thus no long- 
ranger order may exist unless K = 0, which is however 
the limit of infinite density. 

Dynamical response function and Luttinger-liquid 
analysis.- These results can be analyzed within the 
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Luttingcr-liquid theory. We want to calculate the 
imaginary-time r correlation function F(x, r) = 



{T T e 



i2<h(x,T) 



-2i<t>(0,0) 



), on a finite size system of length 



L. It is known [27j from bosonization that F(x, r) = 

(ira/L) 2K /[smh 2 (^) + sin 2 )]* valid in the long- 
time (low-energy) limit ut 3> a where the a is a short- 
distance cutoff of the order of n~ l , it is the velocity of the 
excitations and K the Luttinger exponent. After Fourier 
transforming F(x, r) in q-space with q = 2irj/L, we get: 

(l + e- , j 

[T(j + K)][T(K)T(j + 1)}- 1 2 F 1 (K + j, K;j + 1; e" 2 ^ 

where 2-F1 and T arc the Hypcrgcomctric and Eulcr func- 
tions. 

Using an Ansatz e" TU)li - q ' ] to fit the long-time behavior 
F(q,r -> 00) ~ e - 2 r r (^+J) we get w(q) = 2nuK/L + 
u\q\, where we have used the even parity of the response 
function. Thus, there should be no roton gap at q = 2n 
in the infinite size limit. For finite size, an apparent ro- 
ton gap (vanishing as l/L) can be seen. This gap can be 
traced to the zero mode contribution to the correlation 
functions. All the fits to the RQMC data presented in 



Figs. [2][3] reproduce remarkably well this l/L scaling [32| . 
and are consistent with our previous findings on the den- 
sity dependence of Luttinger- if exponent K (n) 0]. 
Conclusions.- In conclusion, the analysis of the RQMC 
simulational data neatly lead to two main conclusions, 
namely that there are no roton excitations appearing at 
the first star of the reciprocal lattice and that no long 
range order may exist in the whole crossover from the 
TG gas at low density to the quasi-ordered DDW state 
at high densities. The RQMC data analysis is in remark- 
able agreement with what expected for a super-strongly 
(Correlated Luttinger-liquid state, with the inclusion of 
significant dynamical effects. The realization of ID dipo- 
lar quantum (molecular) gases in the TG to the DDW 
r,egimc is within reach of current experimental efforts Q 
and thus our predictions on the excitation spectrum and 
on the absence of the roton minimum, can be tested in 
future experiments by means of e.g. Bragg spectroscopy 
techniques [28| . 
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